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ABSTRACT 

Context. Recent observations demonstrate that the thin and thick disks of the Galaxy have different chemical abundance trends and evolution 
timescales. 

Aims. The relative abundances of a-elements in the thick Galactic disk are increased relative to the thin disk. Our goal is to investigate the 
cause of such differences in thick and thin disk abundances. 

Methods. We investigate the chemical evolution of the Galactic disk in the framework of the open two-zone model with gas inflow. 
Results. The Galactic abundance trends for a-elements (Mg, Si, O) and Fe are predicted for the thin and thick Galactic disks. 
Conclusions. The star formation histories of the thin and thick disks must have been different and the gas infall must have been more intense 
during the thick disk evolution that the thin disk evolution. 
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1. Introduction 

Recent studies of the kinematics and ages of disk stars of our 
Galaxy have revealed the presence of two distinct populations 
in the Galactic disk, called the thick and the thin disk popula- 
tions. Gilmore & Raid (1983) have shown that the stellar den- 
sity distribution towards the Galactic South Pole requires the 
presence of at least two kinematically distinct disk structures. 
For z< 1 kpc their data are fitted by a single exponent and by 
a second exponent in the range 1 - 5 kpc above Galactic plane. 
The authors considered this as evidence of the presence of an 
old thin disk and a thick disk. The thick disk stars have a scale 
height from 760 to 1450 pc (Gilmore & Raid 1983; Robin et al. 
1996; Reyle & Robin 2001; Chen 1997); the scale height of 
the thin disk is 100 to 340 pc (Gilmore & Raid 1983; Chen 
1997). The velocity dispersion (cr u ,o- v ,cr w ) of the thick and 
thin disk stars is also different; it is (63+6,39+4,39+4) km/s 
and (39±2,20±2,20±1) correspondingly (Soubiran et al. 2003). 

Jimenez et al. (1998) used the Hipparcos color-magnitude 
diagram of field stars and derived a minimum age for the 
Galactic disk of about 8 Gyr. Fuhrmann (1998) found that the 
thin disk stars in his sample are younger than 8-9 Gyrs while 
the thick disk star ages exceed 10 Gyrs. Bernkopf et al. (2001) 
suggested that the maximum age for the thin disk stars is near 
9 Gyrs whereas the thick disk stars have ages from 12 to 14 
Gyrs. Bensby et al. (2003) found that the thick disk stars are 
older on average than the thin disk stars, 1 1 .2 + 4.3 Gyrs and 
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4.9 + 2.8 Gyrs. In Feltzing et al. (2003) they derived slightly 
higher average ages - 12.1 ± 3.8 Gyrs and 6.1 ±2.0 Gyrs for 
thick and thin disks correspondingly. Thus, the thin and thick 
disks were formed at different epochs. 

The chemical characteristics of these two disk subsystems 
of the Galaxy are also noticeably different. A kinematic study 
of the local G dwarf metallicity distribution of our Galaxy pro- 
vides evidence that the abundance distribution below [Fe/H]~ 
-0.4 contains two overlapping distributions, the thin and the 
thick disks (Wyse & Gilmore 1995). The mean metallicity of 
the thick disk population, <[Fe/H]>, covers the range between 
-0.5 and -0.7 dex (Soubiran et al. 2003; Wyse & Gilmore 1995; 
Robin et al. 1996) while the mean value of the thin disk star 
metallicity is -0.17 to -0.25 dex (Soubiran et al. 2003; Wyse & 
Gilmore 1995). 

The thick and thin disk abundance trends partly overlap in 
the range -0.8 < [Fe/H] < -0.4 but they are different in [a/Fe]. 
Gratton et al. (1996) pointed out that changes in [Fe/O] values 
of the disk stars can be explained by a systematic difference 
in the chemical composition of the thin and thick disk stars. 
Later Gratton et al. (2000) studied this interpretation in detail; 
as their [Fe/0]vs[0/H] diagram shows, there are two groups of 
disk stars with [O/H] > -0.5 - the thin disk stars with [Fe/O] > 
-0.25 and the thick disk stars with [Fe/O] < -0.25. 

The thick disk has a higher [O/Fe] ratio than the thin disk 
at sub-solar metallicities. The thick disk also shows signatures 
of chemical enrichment by type la supernovae (Bensby et al. 
2004). Fuhrmann (1998) have demonstrated a clear disk sep- 
aration in [Mg/Fe] in the range -0.6 < [Fe/H] < -0.3 in their 
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[Mg/Fe]vs[Fe/H] diagram - the stars with thick disk kinemat- 
ics have [Mg/Fe]^+0.4 while the thin disk stars have [Mg/Fe] 
value decreasing from +0.2 to 0.0 dex. Bensby et al. (2003) 
have analyzed the spectra of F and G disk dwarfs with metallic - 
ities located in the range -0.8 < [Fe/H] < +0.4 and have found 
that the thin and thick disk abundance trends are clearly sep- 
arated at [Fe/H] < 0. Bensby et al. (2005) confirms their pre- 
vious results - there are distinct and separate abundance trends 
between the thin and thick disks. The thick disk stars are more 
enhanced in their a -element abundances than the thin disk at a 
given [Fe/H] below solar metallicities. 

On the other hand, Chan et al. (2000) have studied the 
chemical composition of 90 F and G disk dwarfs and did not 
find a clear separation of [a/Fe]vs[Fe/H] into thin and thick 
disks. As shown by Prochaska et al. (2000), this situation ap- 
pears because the old high - metallicity thick disk stars with 
r e ff < 5700 K were not been included in the sample of dwarfs 
investigated by Chan et al. (2000). 

Recent investigations showed that abundance trends for Mn 
and Eu are also different for these two disk components. In 
particular, a detailed study of the Mn abundance trend of the 
disk and metal-rich halo stars was carried out by Nissen et al. 
(2000). Prochaska & McWilliam (2000) have revisited these 
data (with the new hyperfine structure of the Mnl line) and 
have shown that thin disk stars with -0.8 < [Fe/H] < -0.2 
have [Mn/H]^ -0.1 while thick disk stars [Fe/H] < -0.6 have 
[Mn/Fe]^ -0.3. Mashonkina & Gehren (2000, 2001) report that 
a step-like change exists in [Eu/Ba] and [Ba/Fe] ratios at the 
the transition of the thick to thin disk. They also found that eu- 
ropium is overabundant relative to barium and iron in the halo 
and thick disk stars. Nearly solar [Eu/Fe], [Ba/Fe] and [Eu/Ba] 
ratios are found for thin disk stars. 

In this paper, we consider the thick and thin disks as two 
separate components of the Galactic disk which have been 
formed at different epochs and have different evolutionary 
timescales. 

Matteucci & Greggio (1986) suggested that the evolution 
of the abundance of oxygen vs iron can be explained by an ori- 
gin of iron and oxygen from different SN types. Such an inter- 
pretation is widely accepted since provides an explanation for 
the observed overabundances of a-elements in the halo. Iron is 
predominantly synthesized by type la SN (intermediate mass 
stars in binary systems), while a elements (in our case O, Mg, 
Si) are synthesized by massive stars. Thus, abundance ratios 
depend on the stellar nucleosynthesis, stellar lifetimes and ini- 
tial mass function. But the variation of abundance ratios as a 
function of metallicity or time depends on the star formation 
history as well. This allows us to use the variations of [a/Fe] 
as an indicator of the star formation history (Matteucci 1992). 
Since the [a/Fe] ratios for the thick and thin disk stars are dis- 
tinct, we make the assumption that star formation histories of 
the two disk components must be different. The aim of this 
paper is to find such parameters for the star formation history 
which allow us to reproduce the observed a - element abun- 
dances in the thin and thick disks of the Galaxy at metallicities 
below the solar value. 

This paper is organized as follows. In Sect.2 we describe 
the model and its main components for the chemical evolution 



of the Galactic thin and thick disks. In Sect. 3 we show and 
discuss the model predictions for the abundances of elements 
considered (Mg, Si, Fe) for the thick and thin Galactic disks. In 
Sect. 4 we summarize the main conclusions. 

2. The model 

Our model for the chemical evolution of the Galaxy is de- 
scribed in detail by Pilyugin & Edmunds (1996) and Pilyugin 
(1993). We construct an open two-zone chemical evolution 
model of the Galactic disk in order to investigate the distinc- 
tion in chemical characteristics of the Galactic disk subsys- 
tems. The Galactic disk is divided into two zones - the thin and 
thick disks; the mass of gas, heavy elements, stars and stellar 
remnants are computed as a function of time for each compo- 
nent. The gas infall from the extrahalo takes place during their 
evolution. 

In our model the continuous star formation process in the 
galactic disk is considered as a sequence of bursts with a gen- 
eration of stars formed during each burst. 1 Each stellar gener- 
ation contributes to the chemical enrichment of the interstellar 
medium; the next stellar generation is formed from the galac- 
tic gas enriched by heavy elements ejected from the previous 
stellar generation. We set an interval between bursts of 50 Myr 
as such a value is most suitable for a Galactic disk model. The 
number of bursts depends on the interval between bursts and 
the duration of disk evolution. Thus, during galactic evolution 
the heavy element fraction and number of stars increases while 
the galactic gas mass decreases. We have used use the equations 
(3)-(5) from Pilyugin (1993) for a description of the temporal 
variations of the mass of gas, stars and heavy elements during 
galactic evolution. It is assumed that the evolution time of the 
disk of our Galaxy is 13 Gyrs (Cowan et al. 1991) and the burst 
interval has been set to 50 Gyr (Pilyugin & Edmunds 1996) . 

The star formation rate if/(t) in thick and thin disks is de- 
scribed in the following way (Pilyugin & Edmunds 1996): 

ua (t-e-« T >°r, t<T top 

where T top and T s f r are free parameters of the star formation 
rate. 

The star formation history is described by a set of star for- 
mation bursts whose amplitudes are derived as 

M bi = f " if,(t)dt, (1) 

where and are the times of beginning of j- and j-1 star 
formation bursts. The interval between consistent star forma- 
tion bursts is 50 Myr. 

1 This type of description of star formation rate is usually used for 
a numerical model. In our case "star fomation burst" means a number 
of stars born in a short time interval. The advantage of using this type 
of description is that we can obtain any type of star formation history 
by varying the time interval between bursts and the burst amplitude. 
Such a description of the star formation history is rather general and 
allows us to produce a star formation history in the form of bursts as a 
continuous star formation history. 
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Table 1. Parameters of selected models 



Model 


Thin disk 


Thick disk 


T,„„[Gyr] 


1 


1 


T s/r [Gyr] 


8 


5 


T,„/[Gyr] 


5 


7 


ao 


0.06 


0.1 



0.01 



0.008 - 
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Fig. 1. The star formation and the accretion rates as a function of time 
in the model for thin (solid line) and thick (dashed line) disks. 



We also assume that infall of intergalactic gas takes place 
on the disk during the galaxy's lifetime. According to Pilyugin 
& Edmunds (1996), the infall rate the is described by the func- 
tion 

A(t) = a e-' !T '"f, 

where Ti„f and oq are free model parameters. The infalling gas 
has the primordial chemical composition since, according to 
Tosi (1988), the accretion of a gas with metallicity up to 0.1 
solar abundance gives the same result as accretion of a gas with 
primordial chemical composition. 

The star formation histories of the disk components were 
chosen so that the thick disk stars were formed 10-13 Gyr ago 
while the majority of the thin disk stars would have ages less 
than 10 Gyr. 



1.2 r 




0.001 0.010 0.100 1.000 10.000 100.000 
time, Gyr 

Fig. 2. The cumulative amounts of newly synthesized elements ejected 
by a stellar generation as a function of time. Each curve is normalized 
to its value at t = 13 x 10 9 yrs. 

2.1. Nucleosynthesis 

The abundance of the elements synthesized by a stellar gen- 
eration and ejected to the interstellar medium depends on 1) 
the synthesis of elements by stars of various masses 2) the 
number of stars formed in a given interval of stellar masses 
i.e. on the initial mass function. We adopted the yields of the 
Padova group (Marigo 2000; Portinari et al. 1998)) for calcu- 
lation of the heavy element yields synthesized and ejected by 
massive, intermediate and low-mass stars during and at the end 
of stellar lifetimes. The elements that have been included in the 
calculation are He,C,0,Mg,Si,S,Fe,Ca. The initial mass func- 
tion is described by the Salpeter law with A=2.35 (Salpeter 
1955). In this paper we have used the yields with above-solar 
initial metallicity, Z=0.015. Drastic differences exist between 
predicted yields for stars with initial metallicity Z = and Z > 
(Woosley & Weaver 1995) therefore the effect of a reason- 
able initial metallicity on final stellar yields is negligible for 
stars with Z > 0. 

The Mg yield in the paper of Portinari et al. (1998) was 
determined from Woosley and Weaver's yields (Woosley & 
Weaver 1995) which give an underestimated Mg value (see 
Thomas et al. (2000)); therefore the [Mg/Fe] ratio is not re- 
produced well in their paper (Portinari et al. 1998). 

Francois et al. (2004) tested the Woosley and Weaver yields 
and found that the Mg yield is underestimated. In order to re- 
produce the Mg abundance there is a need to increase the pre- 
dicted Mg yield from stars with masses 1 1 - 20 M Q by a factor 
of 7 while the yields from stars with masses of more then 20 
M Q need be lowered about two-fold. 

Making use of Portinari's data in the calculation of the 
model of chemical evolution we found that stars formed from a 
gas with the initial metallicity Z < 0.15 give a lower Mg yield 
so as to reproduce the observed data. But stars with Z=0.15 
give even higher Mg yields; we had to lower the predicted Mg 
yield from the stars with 9-15 M Q by a factor of 1.5 so that the 
model results would be in good agreement with the observation 
data in the [Mg/Fe] vs[Fe/H] and [Si/Fe] vs[Mg/Fe] diagrams. 

Fig. 2 demonstrates the cumulative mass of magnesium, sil- 
icon, oxygen and iron synthesized by a single stellar genera- 
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tion. The iron yield contributed by SNIa is taken into account 
in the calculation of the total Fe yield. It was assumed that type 
la supernovae evolve according to a SD (single degenerate) sce- 
nario, in the framework of which the accretion of He or H on 
a white dwarf in a binary system is considered (see Nomoto 
1982; Hachisu et al. 1996). The Fe contribution from SNII and 
SNIa types at all times of evolution is about 1/3 and 2/3 corre- 
spondingly. 

3. The model results 

3.1. The a -element abundances 

In this paper we compared the model predictions with the 
observations of abundances of the thin and thick disk stars 
from the paper of Mishenina et al. (2004). Data of Bensby 
et al. (2003) (thick disk star ages) were used in addition for a 
more extended sample of thick disk stars. In our paper, ages 
for Mishenina's data set were calculated using the Bertelli 
isochrones (Bertelli et al. 1994). 

The model predictions for the thin disk relative abundances 
are shown in Fig. 3,5. The best fit model parameters are listed 
in Table 1 . The star formation and accretion rates for the thin 
disk are shown in Fig. 1 a,b. 

The parameters were chosen so that the use of such param- 
eters for the star formation and accretion rate would reproduce 
the observed data. 

As Fig. 3 demonstrates, the use of the above-mentioned pa- 
rameter values allows us to reproduce the [Mg/Fe], [Si/Fe] and 
[Mg/Si] thin disk ratios quite well. The model track agrees 
closely with the line drawn by the least-squares method through 
a cloud of points marking the positions of the thin disk stars in 
the [Mg/Fe] vs[Fe/H] and [Si/Fe]vs[Fe/H] diagrams. However, 
the model poorly reproduces the observed data at the super- 
solar metallicities especially for Si. 

Model predictions for the thick disk are shown in Fig. 4. 
The best fit parameters of thick disk evolution are presented 
in Table 1 . The star formation and accretion rates for the thick 
disk are shown in Fig.l a,b. As Fig.4 shows, the solid line of 
the model prediction for the thick disk is in agreement with 
dashed line obtained by averaging the least-squares method of 
the observed relative abundances of Mg, Si and Fe. 

The model predictions for thin disk element abundances as 
functions of time are worse than for relative abundances. Fig. 5 
shows that the predicted abundance trends for the thin disk are 
slightly lower than Mishenina's observations averaged by the 
least-square method. 

The obtained values of parameters of thin disk evolution 
are in good agreement with parameters of the best fit model 
of the Galactic disk of Pilyugin & Edmunds (1996) who have 
investigated the "age - [Fe/H]" and "age - [O/H]" ratios and 
the solar neighbourhood metallicity distribution function. The 
assumed disk age is 13 Gyr; the obtained age - metallicity re- 
lation reproduces quite well the observed disk star abundances 
of Edvardsson et al. (1993). We compared the model predic- 
tion for the thin disk with the observed abundances of disk stars 
obtained by Edvardsson et al. (1993) for the whole disk. Fig. 6 
demonstrates that under the same parameters the thin disk mod- 
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Fig. 3. The model predictions of the relative abundances of Galactic 
thin disk stars. The data are taken from Mishenina et al. (2004) and 
Bensby et al. (2003)(thick disk star ages only). The asterisks are thick 
disk stars, the small diamonds are thin disk stars, squares are observed 
data averaged in 10 bins with equal numbers of stars. The dashed line 
is a curve drawn by the least-squares method in the observed data. The 
model prediction for the thin disk is indicated by a solid line. 



elled curve is in agreement with the observed age - metallicity 
relation. 

Recently Nordstrom et al. (2004) presented new determi- 
nations of metallicities and ages of F and G-dwarfs in the so- 
lar neighbourhood. We used the probabilities of the stars with 
measured velocities (U, V, W) from the paper of Mishenina et al. 






[Mg/Fe] look-back time, Gyr 

Fig. 4. The model predictions of the relative abundances of Galactic Fig. 5. The model predictions of the abundances of Galactic thin disk 
thick disk stars. Symbols are the same as in Fig.3. stars as a function of time. Symbols are the same as in Fig.3. 



(2004) and the ages and metallicities of Nordstrtim's survey 
for the stars investigated in both papers to build the new "age 
- metallicity" relation, Fig. 7 (top). Fig. 7 shows that predicted 
abundances of the thin disk model (t ga i = 13 Gyrs) are also 
lower on average, as we can see in the case of Mishenina's 
data. 

The computed probability of belonging to thin/thick disk 
for stars with measured velocities (U,V,W) of Mishenina et al. 
(2004) are in good agreement with the ones of Bensby et al. 
(2003). Nordstrom et al. (2004) have compared their photomet- 
ric metallicities with the spectroscopic values for F and G-stars 



from Edvardsson et al. (1993) and found that the agreement is 
excellent with a mean difference of only 0.02 dex and disper- 
sion around the mean of 0.08. 

Thus, there is the question of whether this "overestimated" 
(in comparison with the average metallicity of the whole disk) 
mean value of abundances is usual for thin disk stars. If this is 
the case, it needs to provide a solution to this problem since it 
is not enough to change only star formation and accretion pa- 
rameters in order to reproduce the thin disk age - metallicity 
relation in the framework of our model. Such a age - metallic- 
ity relation (Fig. 5) can be explained either by thin disk pre- 
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Fig. 6. The "age-metallicity" relation of the Galactic disk. Data are 
taken from Edvardsson et al. (1993). The model prediction of the thin 
disk is indicated by the solid line. 

enrichment or an infall of metal-enriched gas. We consider that 
in our model the thin disk is pre-enriched as we have used stel- 
lar yields with above-solar initial metallicity. Tosi (1988) calcu- 
lated the model with metal-rich gas infall for the Galactic disk. 
Table 2 in her paper demonstrates that the effect of a metal-rich 
infall on the age - metallicity relation is negligible in the solar 
neighbourhood. We ran a set of models with infall of enriched 
gas with Z from 0,lZ o to 1 Z Q . Solar proportions for the abun- 
dances of the various metals were assumed. We found out that 
the infall of enriched gas does not noticeably change the posi- 
tion of the modelled curve on the age - metallicity diagram. 

Unfortunately, the number of stars belonging kinematically 
to the thin disk exceeds the number of stars belonging to the 
thick disk; among the latter there are not enough objects for 
which one can reliably determine the ages. Therefore, the trend 
in the "age - metallicity" diagram for the thick disk stars is less 
reliable than for the thin disk stars. The few thick disk stars with 
determined ages were taken from Bensby et al. (2003, 2004) to 
increase the number of thick disk stars in the "age- metallicity", 
"age- [Mg/H]" and " age - [Si/H]" diagrams. 

The situation is slightly complicated by the presence of 
stars with thick disk kinematics and thin disk metallicity in 
the observed data of Mishenina et al. (2004). These stars are 
marked with asteriscs as belonging to the thick disk population 
but are located in the region occupied by a younger and more 
metal-rich thin disk population (e.g. Fig. 3). 

The "age- metallicity" diagram (see Fig.4) shows that the 
thick disk stars are less metal-rich and older than the thin disk 
stars. We have a problem determining the age range of thick 
disk stars. The use of various sets of isochrones gives ages of 
the oldest thick disk stars exceeding the age of Univers by sev- 
eral Gyrs. 

Until the ages of the oldest stars are determined accu- 
rately, we can consider the "age - metallicity" relation of thick 
disk stars as only an estimate of predicted abundance ratios. 
Nevertheless, the predicted "age - metallicity" ratio (as "age 
-[Si/H]" ratio) reproduces the average values of the observed 
data for thick disk stars marked by squares with errors in Fig.4. 

The abundance ratio of Mg and Si depends on the yields 
of single stars with different masses and lifetimes used in the 
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Fig. 7. The "age - metallicity" relation of the thin and thick disk of the 
Galaxy obtained by Nordstrom et al. (2004)(top), and by Mishenina 
et al. (2004)(bottom). The thin disk stars are marked by large dia- 
monds, thick disk stars - asterisks. The solid line is the model predic- 
tion for the thin disk, the dashed line is drawn by the least-squares 
method. 



model. As Fig. 3 and 4 show, the Mg and Si abundance ratios 
are in agreement with the observed data of the thin as thick 
disks, except the observed data at super-solar metallicities. The 
behavior of the thin and thick disk abundance trends at super 
solar metallicities requires further study. 

3.2. The metallicity distribution function 

The successful modeling of the chemical evolution should re- 
produce the G-dwarf metallicity distribution function of the 
Galactic disk. The model prediction was compared with ob- 
servations of the metallicity distribution function in the solar 
neighbourhood (Fig. 9). It is known that 94 % of solar neigh- 
bourhood stars belong to the thin disk whereas the remaining 
6 % belong to the thick disk population (Robin et al. 1996). 
Therefore the observed metallicity distribution functions in 
Fig. 9 were compared with the model prediction for the thin 
disk. With the above-mentioned parameters of the star forma- 
tion and gas accretion in the thin disk (see Table 1), the model 
reproduces the observed metallicity distribution functions quite 
well (Fig.9). 
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Fig. 8. The model predictions of the abundances of Galactic thick disk 
stars as a function of time. Symbols are the same as in Fig. 3. 



3.3. Comparison with previous work and discussion 

A model of chemical evolution should account for the origin 
and characteristics of the subsystems of the Galaxy. A set of 
papers dedicated to aspects of galactic chemical evolution ex- 
ists in the literature. Most of these papers deal with the standard 
open models of chemical evolution with inflow of intergalactic 
gas in the Galaxy. 

Chiappini et al. (1997) considered a two-zone model for 
the disk and the halo, the two-infall model with two episodes 
of intergalactic gas infall, in the formation process of the thin 



0.2 - 



0.1 - 




-1 .0 



-0.5 



0.0 



0.5 



[Fe/H] 



Fig. 9. The metallicity distribution function. The dotted line is the ob- 
served distribution obtained by Hou et al. (1998), the dashed line is 
the observed distribution obtained by Jorgensen (2000), the solid line 
is the model prediction for the thin disk. 



disk and the halo-thick disk. Their model predicts the abun- 
dances of 16 chemical elements except Mg and 3 He. The au- 
thors conclude that the comparison between theory and obser- 
vations suggest a lower iron yield from SNIIs and an increased 
Mg yield from SNIIs. In the framework of Chiappini's model, 
Francois et al. (2004) have computed the chemical evolution 
of 12 elements (including O, Mg, Si) and have reproduced the 
results of observations for O, Si and Ca. Their [Mg/Fe] ratio 
is underestimated owing to the low Mg yield from massive 
stars. Chang et al. (1999) introduced the two-component model 
of chemical evolution which assumes that the thin and thick 
disk of the Galaxy have formed during two accretion episodes. 
Their models reproduce the metallicity distribution function, 
age - metallicity relation and [0/Fe]-[Fe/H] ratio for the Solar 
vicinity. Alibes et al. (2001) analyzed the evolution of all sta- 
ble isotopes between H and Zn. Their model also assumes the 
formation of the Galaxy in two main episodes of exponentially 
decreasing infall. The first episode forms the halo and the thick 
disk from infall of primordial extragalactic gas over 1 Gyr. 
The thin disk forms by infall of gas with metallicity 0.12 Zo 
over 7 Gyr. Their model predictions of abundance ratios are in 
agreement with observation except for Mg (the authors used 
the yields of Woosley & Weaver 1995). 

Goswami & Prantzos (2000) studied the evolution of the 
abundances of the intermediate mass elements ( C to Zn ) in 
the Galactic halo and local disk. Their model assumes strong 
outflow in the halo phase and infall in the disk, which allows 
them to reproduce the corresponding metallicity distribution. 
The evolution of a - elements O, Si, S and Ca is well under- 
stood with the assumption that SNIa contribute most of the 
Fe in the disk but the Mg is underproduced in the Woosley 
& Weaver (1995) yields, therefore their Mg yield is relative to 
the observation in the disk. The model predictions cover a wide 
[Fe/H] range, from -4.0 up to 0.0, but do not show the clear lim- 
its of the model components. One can define the limits only by 
comparing their predictions with observations. Unfortunately, 
in the above-mentioned papers model predictions have not be 
compared with observations for the thin and thick disks sepa- 



8 



T.V. Nykytyuk et al.: The Galactic thick and thin disks: differences in evolution 



rately and we can only conclude about the evolution of the disk 
as whole. 

In contrast, we consider the thick disk evolution separately 
from the halo evolution. The new observations demonstrate the 
clear separation of Galaxy disk stars into the thin and thick disk 
population in the range -0.6 < [Fe/H] < -0.3 (see Introduction) 
that shows the need to consider the Galactic thick disk evolu- 
tion separately, not only from the thin disk evolution but also 
from the halo evolution. 

Ferrini et al. (1992) have considered the evolution of the so- 
lar neigbourhood as subdivided into halo and disk regions. The 
chemical abundances of 14 species (in particular O, Mg, Si, Fe) 
have been predicted for the disk and halo. Their [0/Fe]vs[Fe/H] 
diagram demonstrates the predicted model curve covering the 
[Fe/H] range of -2.5 to 0.5 (0.0 < [O/Fe] < 0.5) for the disk 
and -2.5 to -1. (0.2 < [O/Fe] < 0.5) for the halo. Travaglio et al. 
(1999) have calculated the Galactic chemical evolution of the 
elements from Ba to Eu in the framework of the three-zone 
model of chemical evolution (the halo, the thin disk, the thick 
disk). Later, the Galactic evolution of Pb, Cu and Zn was com- 
puted in framework of the same model (Travaglio et al. 2001; 
Mishenina et al. 2002). As an example of the model results, 
Travaglio et al. (1999) computed the values of [O/Fe] vs [Fe/H] 
in the three Galactic zones of their standard model. We can 
see from their Fig. 3b that the thick disk phase covers the in- 
terval -2.5 < [Fe/H] < -1, and the thin disk phase begins at 
[Fe/H] > -1.5. However, the new observations show clearly that 
the thick disk star population continues up to [Fe/H] =-0.3 at 
least (Mishenina et al. 2004). According to Mishenina et al. 
(2004), star formation in the thick disk stopped when the en- 
richment was [Fe/H]=-0.3, [Mg/Fe]=+0.2 and [Si/Fe]=+0.17. 
Our model predicts that the thick disk relative abundances 
(Fig. 4) come abruptly to an end at [Fe/H] — 0.3 and deviate sig- 
nificantly from the least-square method curve at [Mg/Fe] — 1-0.2 
and [Si/Fe]=+0.17. In our case the thick disk star formation 
stopped at similar values as noticed by Mishenina et al. (2004). 
But the observations show that there are a few thick disk stars 
with [Fe/H] > -0.3 (see also Bensby et al. 2005). Our model 
does not predict the formation of these stars in the thick disk. It 
is possible that these stars were not formed in the thick disk but 
were born in a galaxy - satellite and accreted into the Galactic 
disk. 

Looking at the observed ratios [0/Fe]vs[Fe/H] and 
[Fe/0]vs[0/H] in the thick and thin disks, Bensby et al. (2004) 
concluded that the thin and thick disks were formed in different 
epochs and from a homogeneous gas. In their view, the star for- 
mation rate in the thick disk needs to be rapid to contribute the 
oxygen by SNII at the high metallicities before enrichment by 
SNIa (which do not contribute oxygen) will begin. Yet, the thin 
disk evolution must be slow and the star formation rate low. 
This conclusion is in good agreement with results of the two- 
infall model of Chiappini et al. (1997). Their model predicts a 
timescale of 1 Gyr for the thick disk and halo formation and 
about 8 Gyr for the thin disk formation. Their model success- 
fully reproduced [a/Fe] vs [Fe/H] for the Galactic thick disk- 
halo and the thin disk. But they did not compare their results 
with observed thin and thick abundance disk trends separately 
since at that time such observations had not been obtained yet. 
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Fig. 10. The gas fraction as a function of time in the thin (solid line) 
and thick (dashed line) disk computed by the models. 

An indicator of star formation history in the model is the 
"gas fraction - time" relation since it indicates what fraction of 
galactic gas was converted into stars and during what time it oc- 
cured. We found that the star formation history for the thick and 
thin disks must be different to reproduce the observed abun- 
dance trends. The obtained "gas fraction - time" relation indi- 
cates that the star formation history of the thick disk was rather 
brief as the majority of the thick disk stars was formed 2-3 
Gyrs and had almost stopped more then 10 Gyrs ago (Fig. 10). 
On the contrary, the thin disk stars only began to form 9-10 
Gyrs ago (Fig. 10) and the star formation in the thin disk has 
decreased and is now almost stopped. 

The accretion of intergalactic gas significantly affects the 
abundance predictions for the disk subsystems of our Galaxy. 
The observations show that the thin disk and thick disk exhibit 
parallel slopes of [ar/Fe] vs Fe/H] in the range -0.8 < [Fe/H] 
< -0.3. It is supposed that the slope reflects the contribution of 
different supernovae to the interstellar medium enrichment. In 
our paper, the Fe yield from SNIa is taken into account both 
for the thin and thick disks since the thick disk stars show ev- 
idence of chemical enrichment from SNIa (Mishenina et al. 
2004; Bensby et al. 2003). The difference between the thick and 
thin disk abundance trends is caused by an increased inflow of 
gas mass in our model. The trend of the relative abundances of 
the thin disk stars will change location and approach the region 
occupied by the thick disk stars in the range of -1.0 <[Fe/H] < 
-0.3 if the mass of infalling gas at each accretion episode in- 
creases significantly (more than a factor of 10). In other words, 
the a - element enhancement of the thick disk as compared to 
the thin disk can be caused by more mass inflow in a unit of 
time in the thick disk. 

4. Conclusions 

The chemical evolution of the disk of the Galaxy was investi- 
gated in the framework of the two-zone open model with gas in- 
flow. It was supposed that the Galactic disk is divided into two 
zones - thin and thick disks, differing from each other chemi- 
cally and having different evolution timescales. 

The Galactic evolution of a-elements (Mg and Si) and Fe 
was predicted for the thin and thick disks and was compared 
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with recent observations. Our model predicts that the thick disk 
star formation stopped at [Fe/H]~-0.3 and at [Mg/Fe]~+0.2 
and [Si/Fe]=+0.17; for the thin disk relative abundances our 
model predictions cover the interval -1.0 < [Fe/H] < 0.0, that 
agrees with observations. 

We have obtained the values of parameters of the star for- 
mation history for the thin and thick Galactic disks and con- 
clude that it is necessary to use different star formation his- 
tories for the thick and thin disks to reproduce the observed 
abundance trends. 

Gas infall plays an important role in the appearance of a 
- element enhancement of the thick disk compared to the thin 
disk - the inflow rate must be more intensive per unit of time 
for the thick disk than the thin disk of the Galaxy. 

Acknowledgements. We thank Dr. L.S. Pilyugin, Dr. C. Travaglio and 
especially the referee, Dr. Francois, for helpful discussions and useful 
comments. 
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